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Ole R0mer's method still on the stage. The study of two 
bound eclipsing binaries in quintuple system V994 Her. 
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ABSTRACT 

More than three hundred years ago, Ole R0nier measured the speed of hght only by 
observing the periodic shifting of the observed eclipse arrival times of Jupiter's moons 
arising from the varying Earth- Jupiter distance. The same method of measuring the 
periodic modulation of delays is still used in astrophysics. The ideal laboratories for 
this effect are eclipsing binaries. The unique system V994 Her consists of two eclipsing 
binaries orbiting each other. However, until now it was not certain whether these are 
gravitationally bound and what their orbital period is. We show that the system is in 
fact quintuple and the two eclipsing binaries are orbiting each other with period about 
6.3 years. This analysis was made only from studying the periodic modulation of the 
two periods, when during the periastron passage one binary has an apparently shorter 
period, while the other one longer, exactly as required by a theory. Additionally, it was 
found that both inner eclipsing pairs orbit with slightly eccentric orbits undergoing a 
slow apsidal motions with a period of the order of centuries. 
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1 INTRODUCTION 

Eclipsing binaries are ideal astrophysical laboratories, and 
even more than a century of their intensive photometric 
and spectroscopic monitoring, they still represent the best 
method to determine the masses, radii and luminosities of 
stars. Thanks to modern ground- (and space-) based tele- 
scopes we are able to discover these objects in other galax- 
ies and to apply the same methods as used in our Solar 
neighborhood. One very specific method is the analysis of 
their orbital periods. Using the precise times of minima (cen- 
ters of the eclipses of the components), we can determine if 
the system's period is constant, accelerating, decelerating, 
or periodically alternating. If we detect a periodic shifting 
of the times of minima, we determine that there is an ad- 
ditional component in the system, which is orbiting around 
the barycenter with the eclipsing pair. If the system is mod- 
erately inclined to the observer, the light from the eclipsing 
binary needs more and less time to reach us, as it moves away 
and closer to the observer as it orbits the unseen component. 
This method is in fact the same as used in the 17th century 
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by Ole R0mer when measuring the finite spee d of light usin g 
the eclipse times of Jupiter's moons, see e.g. ICohenI (|l940l ). 



2 THE SYSTEM V994 HER 

Dealing with the eclipsing binaries, we have a few advan- 
tages. First of all, there are currently thousands of eclipsing 
binaries known. Moreover, the time baseline of observations 
for some eclipsing binaries is more than a century long. Im- 
portantly, the observations are very easy to obtain, even 
with small telescopes. 

This is the case of one very interesting eclipsing system 
V994 Her (HD 170314, ADS 11373 AB; V = 7.00 mag). In 
2008 it was discovered (|Lee et all [20081 ) that V994 Her is 
the first (at that time) system consisting of two eclipsing 
binaries. From one point on the sky we can see two different 
sets of eclipses, one with the period of Pa = 2.08 days, while 
the other one has a period of Pb = 1-42 da ys. The star was 
also observed with the Hipparcos satellite (|Perrvman et al.l 
[1997), however the strange eclipsing behavior was missed for 
about 15 years . The complete study of this interesting sys- 
tem was made (|Lee et al.ll2008l ) also on the basis of their new 
spectroscopic observations, yielding a set of physical param- 
eters of all four eclipsing components. Their study showed 



2 P.Zasche and R. Uhlaf 



that the system consists of two pairs: A (B8V + AOV), and 
B (A2V + A4V). AU components are weU-detached and stiU 
located on the main sequence. Both orbits are shghtly ec- 
centric. 

On the other hand, one important question arises, 
whether the two echpsing components comprise one grav- 
itationaUy bounded system, or the system is only an opti- 
cal binary. The mutual orbital period of the two pairs can 
be very long and long-time monitoring is rather time con- 
suming. The system V994 Her also contains one more dis- 
tant component observed interferometrically (|Mason et al.l 
whose period was estima ted of about a fe w thousand 
years. Therefore, the authors (|Lee et al ] I2008') speculated 
that the two echpsing pairs could be identified with these 
two visual components. For a brief review of quadruple sys - 
tems with two eclipsing binaries see lCagas fc Peichal (|2012l ). 
Currently we know o nly six such systems, B V -I- BW Dra, 
¥994 Her, CzeV343 ((Cagas fc Peichal I2OI2I), KIC 424779 1 
l|Lehmann et al.ll2012l ), TYC 3807-759-1 llLohr et al.|[20]l ). 
and OGLE LMC-ECL- 16549 l|Graczvk et al.ll201ll ). see Ta- 
ble □ 



3 ANALYSIS 

Here we introduce an original approach of delay-time varia- 
tions of both eclipsing pairs, showing that the system is in 
fact quintuple and the two inner pairs are orbiting around 
each other with much shorter period. We obtained many new 
observations of both pairs, as wel l as re-analyzed the Hippar- 
cos and AS AS (|Poimanskill2002l ) photometry. The individ- 
ual times of minima are presented in Table lAll These data 
were analyzed simultaneously in a combined approach of fit- 
ting both orbits together, using a well-known method usu- 
ally called "Light-time effect", or "Light-tr avel-time effe ct" 
(hereafter LT TE), described elsewhere, e.g. Ilrwij l| 19591 ) or 
iMaveil (|l990l ). 

This method has been used for dozens of binaries in 
the past, however, V994 Her is the first eclipsing binary 
system, where the method can be applied to both binaries. 
The main advantage of this approach is that both eclipsing 
pairs serve as strictly periodic "stellar lighthouses", whose 
apparent period changes can be studied. 

We used our new code for computing this combined 
approach of deriving both inner orbits (their periods, and 
apsidal motion), together with the long orbit of mutual mo- 
tion of the two pairs. Altogether 15 parameters were fitted, 
using all available times-of-minima observations (for A and 
B pairs: (25+36) published minima together with (18-)-17) 
our new unpublished data points). 

In Figure 1 we plot the O-C diagrams for A and B pairs, 
showing their period changes (the y-axis) with respect to the 
time (the x-axis). In these plots the positive y values indi- 
cate that the detected signal occurs later, while the negative 
values indicate earlier detection than predicted from strictly 
periodic linear behavior. As one can see, the rapid period 
changes near the periastron passage are clearly visible for 
both pairs. Most of the parameters for the LTTE fits for A 
and B are mostly the same. The exceptions are the omega 
angles (the argument of periastron, uja = ljb + 180°), and 
the (O — C) amplitudes {Aa and Ab, see below). The long- 
term modulation arises from the apsidal motion of the inner 
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Figure 1. Plot of the O-C diagrams of both pairs. The dots stand 
for primary, while the open circles for the secondary minima, the 
bigger the symbol, the higher the weight. The red dash-dotted 
lines indicate the LTTE fit, while the black and blue curves rep- 
resent the final fit (LTTE plus the apsidal motion). 




Figure 2. Schematic sketch of the V994 Her system is given (not 
to scale). The orbital planes of both A and B pairs are almost 
perpendicular to the celestial plane (i.e. edge-on to the observer, 
84° and 86°, respectively). 



orbits, because both are slightly eccentric. For the final pa- 
rameters of the fit see Table [2] 



4 RESULTS 

The main result of our analysis is the discovery that the 
two eclipsing pairs orbit around each other and show also 
detectable period modulation together with a slow apsidal 
motion. The period of the motion of apsides for pair A re- 
sulted in about (627 ± 439) years, while for pair B the period 
is about (113 ± 10) years. 

From the parameters of LTTE we determine the semi- 
major axis of the LTTE o rbit, and using the distance 
from the Hipparcos satellite (|Perrvman et al.lll997l ) of tt = 
(3.90 ± 0.74) mas, we also derive the angular separation of 
the two eclipsing binaries on the sky. This resulted in an 
angular separation of Aqi2 = (27.6 ±6.8) mas, which is well 
within the limits for modern stellar interferometers, hence 
its discovery is expect soon. Moreover, this indicates that 
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Table 1. List of currently known double eclipsing systems. 



System Name 


Other designation 


RA J2000.0 


DE J2000.0 


V[mag] 


Period A [d] 


Period B [d] 


Type 


OGLE LMC-ECL-16549 


OGLE LMC-SC3 179761 


05 28 09.41 


-69 45 28.60 


18.216 


164.789640 


0.818033 


EA + EW 


CzeV343 


GSC 02405-01886 


05 48 24.01 


+30 57 03.60 


13.679 


1.209373 


0.806931 


EA + EA 


TYC 3807-759-1 


GSC 03807-00759 


09 30 10.75 


+53 38 59.80 


9.538 


1.305545 


0.227715 


EA + EW 


BV Dra 


HIP 74370 


15 11 50.36 


+61 51 25.25 


8.040 


0.350067 




EW + 


BW Dra 


HIP 74368 


15 11 50.09 


+61 51 41.16 


8.834 




0.292165 


+ EW 


V994 Her 


GSC 02110-01170 


18 27 45.89 


+24 41 50.66 


7.001 


2.083269 


1.420038 


EA + EA 


KIC 4247791 


TYC 3124-1500-1 


19 08 39.57 


+39 22 36.96 


11.645 


4.100871 


4.049732 


EA + EA 



Table 2. V994 Her: The final parameters of the fits for A and B pairs. 







Pair A 




Pair B 




Parameter 


Unit 


Value 


Error 


Value 


Error 


JDo 


HJD 


2448501.1302 


0.0110 


2455375.4555 


0.0062 


P 


Day 


2.0832691 


0.0000038 


1.4200381 


0.0000051 


e 




0.0311 


0.0074 


0.1258 


0.0032 




Deg 


16.0 


4.2 


313.8 


2.8 


duj/dt 


Deg/Cycle 


0.0032 


0.0013 


0.0124 


0.0010 


P3 


Year 


6.33 


0.56 


6.33 


0.56 


To 


HJD 


2456067 


199 


2456067 


199 


A 


Day 


0.0102 


0.0042 


0.0139 


0.0057 


U>3 


Deg 


256.2 


24.1 


76.2 


24.1 


es 




0.747 


0.182 


0.747 


0.182 



the previously mentioned star at ~ 1 arcsec a way is an- 
other star, not the eclipsing one as suggested by iLee et al] 
lj2008h . Membership of this distant component to the eclips- 
ing pairs was suggeste d via similar prope r motions, see e.g. 
the WDS catalogue^, iMason et all l200ll . Thus we have a 
quint uple star system. We currently know only 20 quintu- 
ples l|Eggleton fc Tokovininl i2008l. 

From the parameters of LTTE one can also calculate 
the mass function of the distant body 



/(ma) = 



(ma sin i) 



(mi + m2 + ms) 



Pi 



173.15 A 



see e.g. iMaveij (Il990ll. U sing t he masses of both pairs as 
determined by iLee et al.l (|2008l ). we can also calculate the 
inclination between the eclipsing binary and the LTTE orbit. 
If we label the inclination between the orbit of pair A and 
the LTTE orbit as "m", and vice versa for B, then one can 
derive: 

Aa ms • Al 



sin lA 



mB 
sin i_B 



From this equation and the mass function found from the 
LTTE, we can directly determine the inclination between 
the orbits. This results in ia ~ 37.1° ± 7.3°, while is ~ 
36.8° ±6.8°. Evidently, the inclinations, as derived from both 
A and B pairs, are comparable, and hence we know the 
absolute orientation of the orbit in space. This is the first 
time that the mutual inclination between the orbits of the 
eclipsing pairs has been measured. Only about twenty other 



systems with unambiguous mutual inclination s between the 
eclipsi ng and outer orbits are known, see e.g. iQ'Brien et al.l 
(|201ll ). 



5 DISCUSSION AND CONCLUSIONS 

V994 Her is an interesting target for a future study. How- 
ever, there still remain some open questions. For example, 
these include a detection of the distant component in the 
spectra and determination of its orbital period, etc. Addi- 
tionally, some of the orbital elements of the 6.3-yr orbit are 
still unknown - for instance the longitude of the ascending 
node Q,. The long-term evolution of outer and inner orbits 
should be studied over longer time scales. However, detect- 
ing the mutual motion of the two eclipsing pairs is a unique 
discovery and we can still hope to find similar configurations 
also in other multiple systems. 

Surprisingly, for this analysis there was no need of spec- 
troscopic observations for the radial velocities of the long- 
period system, that would be rather complicated with cur- 
rent observing time allocations on larger telescopes needed 
for such studies. It is noteworthy that all of our new obser- 
vations were carried out with 20-cm aperture or smaller tele- 
scopes by an amateur astronomer. As clearly demonstrated 
by this study, scientifically valuable results can be secured 
with small telescopes by amateur observers. 
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Table Al. Heliocentric minima of V994 Her. 



HJD - 2400000 


Error 


Pair 


Type 


Filter 


48480.29371 


0.00125 


A 


P 


Hp 


48481. 36SS5 


0.00169 


A 


s 


Hp 


52488.4984 


0.0006 


A 


P 


R 


52488.4994 


0.0001 


A 


P 


V 


52488.4997 


0.0003 


A 


P 


B 


52748.9041 


0.0010 


A 


P 


V 


52836.4002 


0.0002 


A 


P 


R 


52937.4701 


0.0001 


A 


B 


R 


53319.71928 


0.00239 


A 


P 


V 


53320.80295 


0.00542 


A 


S 


V 


53904.09882 


0.00169 


A 




V 


53909.26748 


0.00089 


A 


p 


B 


53909.26932 


0.00072 


A 


p 


V 


54290.617 


0.0010 


A 


p 


V 


54314.5204 


0.0008 


A 


s 


V 


54315.522 


0.0010 


A 


p 


V 


54360.3524 


0.0004 


A 


s 


V 


54361.3670 


0.0002 


A 


p 


V 


54383.2729 


0.0003 


A 


s 


V 


54610.3494 


0.0008 


A 


s 


BVR 


54654.08476 


0.00256 


A 


s 


V 


54713.4349 


0.0014 


A 


p 


BVR 


55313.4026 


0.00068 


A 


p 


R 


55314.4827 


0.00046 


A 




I 


55315.4896 


0.00042 


A 


p 


I 


55389.48381 


0.00225 


A 




I 


55413.40460 


0.0010 


A 


p 


VRI 


55641.56985 


0.00012 


A 


s 


R 


55642.66748 


0.00042 


A 


p 


R 


55688.39802 


0.00029 


A 


p 


R 


55691.66016 


0.00011 


A 


s 


I 


55740.47826 


0.00135 


A 


p 


BVRI 


55993.62195 


0.00007 


A 




I 


56041.62826 


0.00101 


A 




I 


56064.44435 


0.00049 


A 




C 


56065.46060 


0.00159 


A 


p 


R 


56138.36968 


0.00059 


A 


p 


I 


56162.36396 


0.00040 


A 




R 


56163.36806 


0.00026 


A 


p 


R 


56187.36470 


0.00076 


A 


s 


IC 


56188.37316 


0.00054 


A 


p 


I 


56210.28311 


0.00180 


A 


s 


C 


56211.28403 


0.00107 


A 


p 


c 


56234.20460 


0.00183 


A 


p 


R 


48479.76947 


0.00305 


B 


p 


Hp 


48480.43493 


0.00263 


B 


s 


Hp 


52471.4363 


0.0004 


B 


p 


UBV 


52381.9864 


0.0010 


B 


p 


V 


52822.20440 


0.00185 


B 


p 


V 


53206.3660 


0.002 


B 


s 


VR 


53504.67119 


0.00454 


B 


s 


V 


53870.18880 


0.00255 


B 


p 


B 


53870.18668 


0.00154 


B 


p 


V 


53887.23816 


0.00216 


B 


p 


B 


53887.23764 


0.00107 


B 


p 


V 


53902.20258 


0.00097 


B 


s 


B 


53902.20047 


0.00048 


B 


s 


V 


53904.28105 


0.00170 


B 


p 


B 


53904.27916 


0.00265 


B 


p 


V 


54283.4131 


0.0001 


B 


p 


V 


54290.616 


0.003 


B 


p 


V 


54298.3787 


0.0002 


B 


s 


V 


54300.4467 


0.0004 


B 


p 


V 


54307.6424 


0.0005 


B 


p 


V 


54332.457 


0.001 


B 


s 


V 


54334.623 




B 


p 


V 


54347.3084 


0.0003 


B 


p 


V 


54364.3439 


0.0003 


B 


p 


V 


54374.2829 


0.0007 


B 


p 


V 


54618.6160 


0.0013 


B 


p 


BVR 


54650.6498 


0.0030 


B 


s 


BVR 


54652.62262 


0.00681 


B 


p 


V 


54653.3886 


0.0013 


B 




BVR 


54699.4777 


0.0014 


B 


p 


BVR 


55375.40806 


0.00037 


B 


p 


I 


55392.44365 


0.0009 


B 


p 


RI 


55392.44779 


0.00055 


B 


p 


I 


55397.49810 


0.00032 


B 




I 


55424.47871 


0.00079 


B 




R 


55654.51846 


0.00097 


B 




C 


55681.50633 


0.00036 


B 




R 


55683.66108 


0.00024 


B 


p 


R 


55691.44630 


0.00013 


B 




I 


55779.48636 


0.00125 


B 




R 


55799.36956 


0.00033 


B 




I 


55821.29667 


0.00095 


B 


p 


RV 


56026.60079 


0.00063 


B 




C 


56048.52317 


0.00028 


B 


p 


R 


56058.46317 


0.00020 


B 


p 


I 


56073.46347 


0.00029 


B 




R 


56102.48296 


0.00058 


B 


p 


R 


56127.42060 


0.00062 


B 




R 


56136.56813 


0.00106 


B 


p 


R 


56154.40166 


0.00039 


B 




C 


56156.44467 


0.00159 


B 


p 


C 


56181.37946 


0.00018 


B 




R. 


56230.28075 


0.00063 


B 


p 


C 
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